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Abstract 

We present a new analytical description of the cellular automaton model for single- 
lane traffic. In contrast to previous approaches we do not use the occupation num- 
ber of sites as dynamical variable but rather the distance between consecutive cars. 
Therefore certain longer-ranged correlations are taken into account and even a mean- 
field approach yields non-trivial results. In fact for the model with v max = 1 the exact 
solution is reproduced. For v max = 2 the fundamental diagram shows a good agree- 
ment with results from simulations. 
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Despite a large number of publications about the cellular automaton approach to traf- 
fic flow (see e.g. |l[ and references therein) only a few of those deal with a systematic 
analytical description. Most works make use of large-scale computer simulations 
which can be carried out very efficiently for this class of models. Nevertheless, ana- 
lytical results - exact or approximate - may give important information relevant for 
a complete understanding of those models. 

The most important exact result is certainly the solution of the model for v max = 10. 
This result has been obtained using n-cluster-approximation |2|, H, i.e. an improved 
mean-field theory taking into account correlations between n neighbouring sites. For 
Vmax = 1 already the 2-cluster approximation is exact ||. For higher velocities the 
n-cluster-approximation for small n already yields very good results for the so-called 
fundamental diagram (flow-density relationship) || |3|] . 

In computer simulation studies there are in principle two different approaches || 
called site-oriented and car-orientecf]. In the site-oriented approach the state of the 
system is specified by storing the state of each cell which can either be empty or 
occupied by a single car with velocity v — 0,1, ... , v max . In the car-oriented approach 
on the other hand one stores the velocity of each car and the distance to the next 
car ahead. 

Since in the cluster-approximation corresponds to a site-oriented approach this anal- 
ogy inspired us to investigate an analytical description based on the car-oriented 
approach, the so-called car-oriented mean-field theory (COMF)[]. The COMF al- 
ready takes into account some longer-ranged correlations so that one can hope that 
it yields at least a good approximation. 

For completeness we briefly repeat the definition of the CA model for single-lane 
traffic flow || in the following. The street is divided into L cells of a certain length 
(for realistic applications 7.5 meters) which can be occupied by at most one car 
or be empty. The cars have an internal parameter ('velocity') which can take on 
only integer values v = 0, 1, 2, ... , v max . The dynamics of the model are described 
by the following update rules for the velocities and the motion of cars ||. In the 
first step all cars with velocities Vi < v max are accelerated by one velocity unit, 
Vi — > v[ = Vi + 1. The following step describes the slowing down due to other cars 
and prevents accidents. All cars with velocities v\ > di (where di is the number of 
free cells in front of car %) decelerate to velocity v" = di. The last step in the velocity 
update is a randomization effect taking into accout several aspects of the driver's 
behaviour, e.g. fluctuations in driving style, overreaction at breaking, and retarded 
acceleration: Every car with velocity v" > will slow down one unit with probability 
p, i.e. v" v" 1 = v'l — 1. In the final step the car then moves v"' sites. These four 
rules, referred to as step 1 to step 4 in the following, are applied to all cars at the 
same time (parallel dynamics). 

tin H this approach has been called particle-oriented. 

^ A brief account of some preliminary results has already been given in Q. 
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COMF for v max = 1 

We denote the probability to find at time t (exactly) n empty sites in front of a vehicle 
by P n (t). As in |2], ||] we change the order of the update steps to 2-3-4-1. This change 
has to be taken into account when calculating the flux f(c,p). It has the advantage 
that after step 1 there are no cars with velocity 0, i.e. all cars have velocity 1. The 
time evolution of the probabilities P n (t) can conveniently be expressed through the 
probability g(t) (g(t) = 1 — g(t)) that a car moves (does not move) in the next 
timestep. 

In order to find the time evolution of the P n (t) we first determine from which con- 
figurations at time t a given state at time t + 1 could have been evolved under the 
rules 2-3-4-1. Take for instance a car - called second car in the following - which has 
n > 1 free sites in front, i.e. its distance to the next car ahead (called first car in the 
following) is n + 1 sites. This configuration might have evolved from four different 
configurations at time t, depending on whether (i) both cars moved in the timestep 
t — > t+ 1 (which happens with probability qg(t)), (ii) both cars did not move (with 
probability pg(t)), (Hi) only the first car moved (with probability pg(t)) or (iv) only 
the second car moved (with probability qg(i)). This means that the second car at 
time t had either n free site in front (cases (i) and (ii)), or n — 1 free sites (case (Hi)), 
orn + 1 free sites (case (iv)). 

The special cases n = 0, 1 can be treated in a analogous fashion. In this way one 



obtains the time evolution of the probabilities as 

Po(t + l) = g(t)[Po(t)+QPi(t)\, (1) 

Px(t + 1) = g(t)P (t) + [qg(t)+pg(t)]Pi(t)+(19(t)P2(t), (2) 

P n (t + 1) = pg(t)P n ^(t) + [qg(t)+pg(t)}P n (t) + qg(t)P n+1 (t). (3) 



A car will move in the next timestep if there is at least one empty cell in front of it 
(probability J2 n >i Pn(t)) an d if it does not decelerate in the randomization step 3) 
(probability q — 1 — p). Therefore the probabilities g(t) and g(t) are given by 

g(t) = qJ2Pn(t) = q[l-Po(t)} 

n>l 

g(t) = P (t)+pY,Pn(t)=P + qP (t) (4) 

n>l 

where we have used the normalization 

E PnW = 1- (5) 

ra>0 

The probabilities can also be related to the density c = N/ L of cars. Since each car 
which has the distance n to the next one in front of him 'occupies' n + 1 sites we 
have the following relation: 

5> + i)P B (f) = i. (6) 

n>0 C 
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Here we are mainly interested in the stationary state (t — > oo) with lim^oo P n (t) = 
P n . In order to determine the probabilities in the stationary state we introduce the 
generating function 

oo 

P{ Z ) = Y J PnZ n+l . (7) 
n=0 

After multiplying the corresponding equation in (|])-(|3]) by z n+1 and summing over 
all equations one finds an explicit expression for the generating function, 

P(z) = qKy ^ y 1 °. 8 

qg - pgz 

The normalization condition (||) and the density relation (^j imply that the generat- 
ing function has to satisfy 

P(l) = 1 and P'(l) = -, (9) 

where P'{z) denotes the derivative of P{z). 

Using (d) it is easy to obtain the probabilities explicitly: 



2qc- 1 + yJl-Aqc(l 
Po = 



Pn = k(eiY = \( p^-Pq) V (B>1) (10) 

where we already used ([]) to express P through the density c of cars. 

To obtain the fundamental diagram we have to calculate the flux. It is given by 
f(c,p) = eg = qc(l — Pq) from which one recovers the exact result 0, |3| 



T - Aqc(l - c) 
f(c, P ) = * ^ (11) 

In 0, Q we expressed the exact solution in terms of the pair probabilities P(nj,nj+i) 
to find two neighbouring sites j and j + 1 in the state (nj,rij+i). Here rij = 
denotes an empty site and rij — 1 a site occupied by a car (with velocity In 
@, |] it was shown that probabilities for larger clusters factorize, i.e. P(ni, . . . , Ul) = 
Y\j=i P{ n ji n j+i)- The 2-cluster probabilities are related to the P n through 

P(l,l) = cP , 

P 2 (1,0) = c(l-c)P h (12) 
P(0,0) = (l-c)^i (n>l). 



^Cars with velocity do not exist after the acceleration step. 
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The factors c and 1 — c appear due to the different normalization of the P n and 
P(nj,rij + i). The P n are normalized by the number of cars whereas the P(jij,rij+i) 
are normalized by the number of sites. 

The fact that the COMF yields the exact result comes not unexpected since the 
COMF takes into account all relevent correlations for the case v max = 1 where only 
nearest-neighbour correlations are non-trivial || [|[ . Although the model for v max > 1 
is know to exhibit a different behaviour |3j the COMF might yield good results in 
this case because it takes into longer-ranged correlations. In the following we will 
investigate the case v max = 2. 



COMF for v r 



The be treated in a similar way as v max = 1. However, it is 

now necessary to introduce two different functions P n (t) and B n (t) describing the 
probabilities to find exactly n empty sites in front of a car with velocity 1 and 2, 
respectively. Proceeding as in the case v max = 1 we find the evolution equations. For 
the stationary state the probabilities obey the equations 

P = go[Po + B ], (13) 

P x = g 1 [P + B ]+pg [P 1 +B 1 ), (14) 

P2 = g 2 [Po + B }+ P g 1 [P 1 + B 1 ]+ P g P 2 , (15) 

^3 = pg2[Pi + B 1 ]+pg 1 P 2 +pg P 3 , (16) 

Pn = pg 2 Pn-2 + PgiPn-1 + PgoPn, (™ > 4) (17) 

and 

Bo = qg [Pi + B 1 + B 2 ], (18) 

By = qg 1 [P 1 + B 1 + B 2 )+g [qP 2 +pB 2 + qB 3 ), (19) 
B 2 = qg 2 [Pi + Pi + B 2 ] + 9l [qP 2 + pB 2 + qB 3 ] + g [qP 3 + pB 3 + qB 4 ] , (20) 
B n = g 2 [qP n -i + pP n -i + qB n ] + g x [qP n + pB n + qB n+1 ] 

+9o [qPn+i + pB n+1 + qB n+2 ] . (n > 3). (21) 

The probabilities g a that a car moves a sites (a = 0, 1, 2) in the next timestep are 
given by: 

g = P + B +pJ2Pn + pBi, 

n>l 

91 = q^Pn + qBt+pJ^Pm (22) 

n>l n>2 

92 = qJ2 B n- 

n>2 

We just mention here that it is possible to derive the identities go = J2n>o Pn and 
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9i + 92 = H n >o B n from (|22|) and (|T3]-pT|) . Using the normalization 



J2 [Pn + B n ) = 1 

n>0 

we have #0 + 9i + 92 = 1. 

The conservation of density leads to the constraint 



(23) 



£>+l) [P n + B r , 



(24) 



n>0 



We introduce the generating functions 



Y n+1 



n=0 

oo 



B(z) = J2 B nZ U+1 - 



(25) 
(26) 



n=0 



These functions have to satisfy P(l) + P(l) = 1 and P'(l) + P'(l) = A due to (|| 
and (El), respectively. 



After multiplication with z n+1 and summation over all the equations ( |i~3| - |21~|) one 
finds 



P(z) 
B(z) 



9W 



l-pg{z 

zg{z) 



(qP + B )z + pB lZ * 



(27) 



z 2 - (q + pz)g(z) 



qP{z) - qB - (qP + pB + qB x )z + (q - p)B lZ < 



(2? 



where we have introduced the function g(z) = g + giz + g2Z 2 . Note that this function 
satisfies g(l) = 1 and g'{l) = -f(c,p) is just the average velocity of the vehicles. 

If one expresses the sums appearing in ( P2"|) by P(z = 1) and B(z = 1) one obtains 
the following relations 



Bi = %o--(B + qPo), 
p p 

9l = p(l-P ) + (g- Vo + (l--)£i, 

q q 

g 2 = q(l-P )-(l + q)B -B 1 . 



(29) 

(30) 
(31) 



With these relations also the normalization condition P(l) + P(l) = lis satisfied. 
Now we can express the generating function completely in terms of the two proba- 
bilities Pq and B only, since go = P ^ B from ([13]). 
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At this point it is surprising that we are still left with two unknowns Pq and Bq 
since we only have one free parameter, the density c. In the following we determine 
a relation between Pq and Bq from analytic properties of B(z). Thus the generating 
functions depend only on one free parameter, e.g. Pq, which then can be related to 
the density via P'(l) + B'{1) = \. 

The denominator of B(z) can be rewritten as z 2 — (q + pz)g(z) = pgi{\ — z){z — 

s + )(z — sJ) where s± are given by s± = ^ 

denominator of B(z) has three zeroes, z = 1 and z = s±. Two of these are located 
in the unit circle since |s_| < 1. These zeroes have to be cancelled by corresponding 
zeroes of the numerator since B(z) has to be analytic in the unit circle (otherwise 
one would not have lim^oo B n = 0). It is easy to see that the numerator indeed has 
a zero at z — 1. Demanding that it also has a zero at z = S- we find the missing 
relation between P Q and B , 

gP(s_) - qB - {qP + pB + qB x ) 8 - + (q - p)B lS 2 _ = 0. (32) 

Owing to (HU) we can regard the generating functions P(z) and B(z) as functions 
of the density c only. The fundamental diagram can then be obtained using the 
following expression for the flux: 

f(c,p) = c[ gi + 2g 2 }. (33) 

Results are shown in Fig. 1 for p = 0.1 and Fig. 2 for p = 1/2. For p = 0.1 
we find an excellent agreement of the two curves. For p = 0.5 we still find an 
excellent agreement for small densities (c < 0.2) and high densities (c > 0.5). Only 
near the maximum there are deviations. For comparision we also have included 
in Fig. 2 results from the site-oriented approach, i.e. the n-cluster approximation. 
The COMF result is much better than the 2-cluster result and comparable to the 
3-cluster approximation. It seems that the COMF tends to overestimate the flux 
whereas the n-cluster approximation yields a lower bound for the flux || . The result 
for p = 0.1 shows, however, that the COMF result is not systematically larger than 
the simulation results. 

For small densities the average distance between the cars is large. Therefore corre- 
lations between cars and neighbouring empty sites are much more important than 
those between two cars. These correlations are better described by the COMF which 
is the reason why the COMF is superior to the cluster approach for small cluster 
sizes in this regime. 

The agreement between simulations and the COMF result is very good up to densi- 
ties close to the critical density c crit ~ 0.19 for p = 1/2 0. Here correlations between 
distances become important. It is to be expected that a combination of the cluster 
approach and the COMF will give a much better agreement of the calculated funda- 
mental diagram and the simulations in the region near the critical point and the flux 
maximum. 



1 ± Jl 



^a 2 

V92 



with a = Thus the 



go+Q9i 
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We also have applied the COMF to generalizations of some modified models |7|, 1. 19| . 
The models discussed in ]7|, || have different modified acceleration rules which for 
Vmax = 1 break the 'particle-hole' symmetry, i.e. the fundamental diagrams are no 
longer symmetric with respect to c = 0.5. It turns out that the COMF gives are 
good agreement with simulations, but is no longer exact (even for v max = 1). A full 
account of these results will be given elsewhere. 
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Figure 1: Fundamental diagram for v max = 2 and p = 0.1. The comparison of the 
COMF result (full line) with results from computer simulations (•) shows an excellent 
agreement. 
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Figure 2: Fundamental diagram for v max = 2 and p — 1/2. The full line is the COMF 
result. For comparison the results from computer simulations (•) and the n-cluster 
approximation for n — 1, 2, 3 (broken lines) are also shown. 
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